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Abstract 

We theoretically discuss the magnetic-field-angle dependence of the zero-energy density of states (ZEDOS) in super- 
conductors. Point-node and line-node superconducting gaps on spherical and cylindrical Fermi surfaces are consid- 
ered. The Doppler-shift (DS) method and the Kramer-Pesch approximation (KPA) are used to calculate the ZEDOS. 
Numerical results show that consequences of the DS method are corrected by the KPA. 
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In the last decade, experimental techniques for the 
magnetic-field-angle dependent measurements of the 
thermal conductivity and specific heat have been devel- 
oped to investigate the superconducting gap anisotropy 
in various materials (jl |2D. The Doppler-shift (DS) 
method [3] has been frequently utilized to analyze such 
. Theoretical development has been 
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cently, we have developed a new method called the 
Kramer-Pesch approximation (KPA) lfl3lll5ll . where the 
contribution of a vortex core neglected in the DS method 
is taken into account 1 13]. Then, the KPA may be a use- 
ful tool that enables a more quantitative analysis. 

In this paper, we show several numerical results for 
the zero-energy density of states (ZEDOS) obtained 
by the DS method and the KPA to demonstrate how 
much the KPA improves quantitative results of the 
DS method. The specific heat over the temperature, 
C(T)/T, is proportional to the ZEDOS in the low tem- 
perature limit T — > 0. 

The expression for the densityjaf states (DOS) N(E) 
in the DS method is given as 
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where 5E = mv F ■ v s is the DS energy, m, v F , and 
v s are the electron mass, the Fermi velocity, and the 
circulating superfluid velocity around a vortex, respec- 
tively. The field-angle dependence is taken into ac- 
count via v s , which is perpendicular to the magnetic 
field H (i.e., v s _L H). Around a single vortex, |v s ] = 
h/lmr (r is the radial distance from the vortex center). 
(• ' ')ds = JC " rdr f Q da - ■■, is the spatial integration 
around the vortex in the cylindrical coordinates (r, a, z) 
with z || H. Here, £o is the coherence length and r a is 
the cutoff length with r a /ft = ^H c2 /H, (H c2 = O /tt^, 
Oo = nr^H, and (Do is the flux quantum). dS F is an 
area element on the Fermi surface (FS) [e.g., dS F = 
fcj: sin 6d<pd6 for a spherical FS in the spherical coordi- 
nates (k, <p, 6), and dS? = kp a t,dd)dk c for a cylindrical FS 
in the cylindrical coordinates (k a b,<p,k c )]. The pair po- 
tential is A s A()A(£pX where Ao is the maximum gap 
amplitude and A(£p) represents the gap anisotropy on 
the FS. A and Vp are functions of the position £p on the 
FS. 

On the other hand, the DOS in the KPA is given 
as |13] 



N(E) = 



v fo t] I CdS F 4 cosh O/6>)] 



I 



|v F | {E-Erf + jf / KPA 

October 21, 2009 



(2) 



Here, A = |A|, and (■ • -)kpa = £" rdr j^" ■ ■ ■ da/nr^ is 
the real-space average around a vortex in the cylindri- 
cal coordinates (r, a, z) with z\\ H. x = rcos(a - 9 V ), 
y = r sin(ff - 9 V ), and E y = AoA 2 y/£;oh. 9 v (k F , a M , 9 M ) is 
the angle of Vf_l in the plane of z = 0, where a and 9 V are 
measured from a common axis lfl2[ IT5I1 . Vf± is the vec- 
tor component of Vp{kp) projected onto the plane normal 
to H — (q-m, 6m)- Here, the azimuthal (polar) angle of H 
is aM (9m) hi a spherical coordinate frame fixed to crys- 
tal axes. |v F x(&f,<*m,0m)I = VFn(gM, 9 M )h(k F , a M 9m) 
and v F o is the FS average of |v F J represented as 012I1 . 
vfo = f dSp\vp±\/ J dSp. The coherence length is ex- 
pressed as £0 = Hv-po/nAo. The smearing factor is set as 
77 = 0.05A . 

We numerically calculate the ZEDOS N(0) for the 
following four combinations: the line-node or the 
point-node gap on the cylindrical or the spherical FS. 
Here, the line-node gap is the af-wave pairing A = 
cos 2cf> sin 2 9 for the spherical FS and is A = cos 2(f> for 
the cylindrical FS. The point-node gap is the s + g wave 
A = (1 + sin 4 9 cos 4-cf))/ 2 for the spherical FS and is 
A = (1 + cos 4 (fc c /2) cos 40) /2 for the cylindrical FS. 

We show the results in Figs. 1 and 2. As shown in 
Fig. 1(a) [the line-node gap on the spherical FS], a cusp- 
like minimum does not appear and the oscillation am- 
plitude is of the order 3.5 % for the KPA, coinciding 
well with a full numerical result in Ref. @|. For the 
DS result, however, the amplitude is of the order 4.9 % 
with less coincidence. As seen in all figures, the DS 
method is liable to give rise to an inflated oscillation 
amplitude. While a cusp-like minimum appears both 
for the DS method and the KPA in Fig. 2(a) [the line- 
node gap on the cylindrical FS], only the DS method 
gives rise to a relatively cusp-like structure in the case 
of the point-node gap on the cylindrical FS as noticed in 
Fig. 2(b), which we mentioned in Ref. lfl3ll . The KPA 
is a perturbation theory, and the DS method corresponds 
to its zeroth order one with a spatially constant gap 01311 . 
Therefore, that discrepancy in Fig. 2(b) signifies an im- 
provement due to the KPA. 

Both the DS method and the KPA are theories valid 
in low-field and low-T region. Taking into account a 



contribution of a vortex core B13I1 . the KPA yields more 
quantitative results with the same computational load 
compared to the DS method. 
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Figure 1 : Azimuthal field-angle ajyi dependence of the ZEDOS for the 
spherical FS. The polar angle 6u = n/2. The cutoff length r„ = 7fn. 
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Figure 2: Azimuthal field-angle cim dependence of the ZEDOS for the 
cylindrical FS. The polar angle 9m = n/2. The cutoff length r a = 7fo- 
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